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Abstract 

In a recent work we introduced a nonlocal version of the Nambu-Jona-Lasinio(NJL) model that 
was designed to generate a quark self-energy in Euclidean space that was similar to that obtained 
in lattice simulations of QCD. In the present work we carry out related calculations in Minkowski 
space, so that we can study the effects of the significant vector and axial-vector interactions that 
appear in extended NJL models and which play an important role in the study of the p, uj and a± 
mesons. We study the modification of the quark self-energy in the presence of matter and find that 
our model reproduces the behavior of the quark condensate predicted by the model-independent 
relation (qq) p = {qq)o(l — o^Pn I f^ m \ + '"')> where &n is the pion-nucleon sigma term and pN 
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is the density of nuclear matter. (Since we do not include a model of confinement, our study is 



Q_i' restricted to the analysis of quark matter. We provide some discussion of the modification of the 

above formula for quark matter.) The inclusion of a quark current mass leads to a second-order 



phase transition for the restoration of chiral symmetry. That restoration is about 80% at twice 
nuclear matter density for the model considered in this work. We also find that the part of the 
quark self-energy that is explicitly dependent upon density has a strong negative Lorentz-scalar 
term and a strong positive Lorentz- vector term, which is analogous to the self-energy found for the 
nucleon in nuclear matter when one makes use of the Dirac equation for the nucleon. In this work 
we calculate the nucleon self-energy in nuclear matter using our model of the quark self-energy and 
obtain satisfactory results in agreement with the values of the scalar and vector nucleon potentials 
in matter found in either theoretical or phenomenological studies. 

PACS numbers: 12.39.Fe, 12.38.Aw, 14.65.Bt 
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I. INTRODUCTION 



In recent years there has been a great deal of interest in understanding the properties 
of quark matter at relatively low temperature and high density [1-5]. Since it is difficult to 
study QCD at finite chemical potential using lattice simulations, the model of choice for such 
studies has been the Nambu-Jona-Lasinio model [6-8]. Of particular interest is the prediction 
of diquark condensates and color superconductivity at high densities. (Such studies may be 
relevant to the properties of neutron stars.) Calculations made for dense matter using the 
NJL model are carried out in Minkowski space. These calculations are limited in that they 
do not include a model of confinement and are, therefore, unable to provide a comprehensive 
description of the hadronic phase present at low density and temperature. However, it is 
generally believed that a good deal of information may be gained by studying quark matter, 
with a proper study of the hadronic phase deferred until some future time. 

As is well known, the NJL model provides a microscopic dynamical description of chiral 
symmetry breaking with the generation of associated quark vacuum condensates and con- 
stituent masses. In the standard versions of the NJL model [6-8], the constituent quark mass 
that is generated in the model is a constant. However, it is known from lattice simulations 
of QCD that the constituent mass goes over to the current quark mass when the quark 
momentum, p 2 , is less than about -2 GeV 2 [9]. It is our belief that, if we are to use the NJL 
model to study dense matter, it is desirable to make the model as realistic as possible. To 
that end, we have introduced a nonlocal version of the NJL model [10] that is able to repro- 
duce the Euclidean-space behavior of the quark mass seen in lattice simulations of QCD [9]. 
To carry out that program we have introduced a momentum-dependent qq interaction in the 
calculation of the quark self-energy and have separated the regularization of the model from 
the specification of that interaction. (This procedure requires the introduction of additional 
parameters into the model.) 

Recently, we have seen an attempt to obtain the quark self-energy in Minkowski space 
by analytic continuation of a Euclidean-space form based upon gluon exchange enhanced 
at small momentum transfers to simulate confinement [11]. The resulting Minkowski-space 
values exhibit resonant-like structures and very large values of the constituent mass. Such 
results do not appear to have any natural physical interpretation. 
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We note that the quark self-energy may be written as 

S(/c 2 ) = A(k 2 ) + B(k 2 )ft 



(1.1) 



in vacuum. In matter there is another four-vector, 77^, that describes the motion of the 
matter rest frame. It is useful to put rf = 1 and to note that, if we work in the matter rest 
frame, we can put rf = [1, 0, 0, 0]. In matter, we have 

Z{k 2 , V -k) = A{k 2 , V -k) + B{k 2 , V -k)J + C{k 2 , V -k)i) , (1.2) 

where we have found it useful to introduce the four-vector 

W = k» - (k "n)rf . (1.3) 

Note that k° = in the matter rest frame. In that frame, we may write 

E(ifc , k) = A(k°, k) - B(k°, fc)7 • jfe + 7°C(A; , k) . (1.4) 

As we will see, A and B satisfy coupled nonlinear equations, while C may be calculated 
independently. (Note that A and Chave the same dimension, while B is dimensionless.) 

We first discuss our results for the case in which we neglect the dependence of A, B, 
and Con k°. We anticipate that the dependence on k° will be relatively weak and justify 
that assumption later in this work. As a further simplification, we will at first neglect B 
and study the behavior of A(k, p), where p is the density of quark matter which is taken to 

— * 

contain equal numbers of up and down quarks. In this case, the maximum value of A(k,p) 

— * 

is found at k — 0, leading to a simpler presentation of our results. We then go on to the 

— * — * 

consideration of the coupled equations for A(k,p) and B(k,p). As usual, we may introduce 
a density and momentum-dependent mass defined by 

M(k,p)= A{Kp ) . (1.5) 
1-B(k,p) 

We provide values of A(k,p), M(k,p) and C(k,p) in the following sections. 

We note that the inclusion of C(k,p) precludes the passage to Euclidean space. The 
analogous problem arises when one introduces a finite chemical potential. As we will see, 

— # 

C(k : p) is quite large at finite density and can only be neglected at very small densities. 
In vacuum, Lorentz invariance leads to a relation between C{k 2 ) and B{k 2 ). However, our 
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formalism does not maintain Lorentz invariance. (For example, our regulator depends only 
upon \k\.) Therefore, in the following we will only calculate the contribution to C(k, p) from 
the matter, which takes the form of two Fermi seas of up and down positive-energy quarks 
with Fermi momentum kp. 

The Lagrangian of our model is 

£ = q(i$> - m°)q + ^ ]T[(gA*g) 2 + (^7 5 A l g) 2 ] 

i=0 

-^E[(^ 7 ^) 2 + (gA l 757,g) 2 ] 

i=0 

+ ^{det[g(l + 75 )g] + det[g(l - 75 )<?]} 

+Aonf (1.6) 

This Lagrangian was used in Ref. [10] with Gy = 0. There, both the condensate values 
and self-energies were given for the up (down) and strange quarks. For the present work 
we neglect the 't Hooft interaction (Go = 0) and our model of confinement. We also 
drop the strange quark from consideration, so that the quark current mass matrix is m° = 
diag(m°,m°) with m" = m° d . Since the 't Hooft interaction contributes to the self-energy 
[10], its neglect leads us to use a larger value of G s than that we would use if the 't Hooft 
interaction were included in our analysis. 

The organization of our work is as follows. In Section II we present the coupled nonlinear 
equations that determine A(k,p) and B(k,p) and also provide an expression for C(k,p), 
and the density- dependent condensate (uu) p = (dd) p . We present values of A(k,p) and the 

— * — * 

condensate, when we neglect B(k,p). In Section III we consider finite values of B(k,p) and 

— » — » — » 

provide results for A(k,p), M(k,p), C(k,p) and the condensate. In Section IV, we consider 
the dependence of A, B and C on k° and provide some results for A(k , k) in vacuum. In 
Section V we use the quark self-energy calculated here to obtain the nucleon self-energy in 
nuclear matter in a simple model. Finally, Section VI contains some further discussion and 
conclusions. 

II. THE QUARK SELF-ENERGY 

In our earlier work we obtained the quark self-energy from the solution of the equation 
depicted in Fig. 1 [10]. There, the open circle is a momentum-dependent quark interaction, 
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FIG. 1: The equation for the quark self-energy that was solved in Euclidean space in Ref. [10] is 
shown. Here, m° is the current quark mass. The open circle represents the momentum-dependent 
qq interaction of the nonlocal model. The third term on the right-hand side of the figure represents 
the 't Hooft interaction and the fourth term arises from our model of confinement. The heavy lines 
are quark propagators which include the self-energy, in their definition. [See Eq. (2.1).] 

obtained by the replacement 

G s ^f(k-k')G s f(k-k') : (2.1) 

where k and k ' are the quark momenta entering (or leaving) the interaction. We have used 

f(k-k') = exp[-(k-k') 2n /2p} (2.2) 

with n = 4 and j3 = 20 GeV 8 in Ref. [10]. The corresponding nonlocal Lagrangian is given 
in Ref. [10]. On the right-hand side of Fig.l, the 't Hooft interaction (third term) and the 
confinement interaction (fourth term) are neglected for the purposes of this work. In the 
second term we have contributions from the negative-energy states in vacuum as well as the 
positive-energy states at finite density, specified by the quark Fermi momentum, kp, of the 
up and down quark Fermi seas. 
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In general, the quark propagator is 



which we will write as 



iS(k°, k) = =, =, =, =, , (2.4) 

(A; - C(k°, k))-f - (1 - B(k°, fc))7 -k - A(k°, k) + ie 

in the matter rest frame. In a first approximation we write 

iS(k°,k) = - % —^ - - , (2.5) 

(A; - C(fc)h - (1 - B(fc))7 -k - A(k) + ie 

— * 

with A(k), B{k) and C(k) density-dependent, in general. (On occasion we will write A(k, p), 
etc.) We see that the presence of C(k) precludes the passage to Euclidean space that was 
made in Ref. [10]. That is analogous to the problem created by the introduction of a chemical 
potential in the formalism. Note that 

(A; - C(k)) 2 - (1 - B{k)fk 2 - A 2 (k) + te 

= [k° -E + {k)+ie][k° -E-{k)-ie\ (2.6) 

with 



E + (k) = C(k) + \]k 2 {l- B{k)) 2 + A 2 (k 2 ) (2.7) 



and 



E~(k) = C(k) - Jk 2 (l - B(k)) 2 + A 2 (k 2 ) . (2.8) 



Here E + (k) may be interpreted as the (on-mass-shell) energy of the positive-energy states, 

— * 

while E~{k) refers to the negative-energy states. 

In order to simplify the notation somewhat, we write 

iS(k°,k) = , (2.9) 

H(A; , k) - A(k) + ie 

where we have defined a four-component quantity 

IP(A; , k) = [k° - C(k), (1 - B(k))k] . (2.10) 
We also introduce a scalar quantity 



V (2*)" 



p s (k) = iN e (-l) jWi^O , (2 .U) 



where the minus sign is due to the closed Fermion loop in Fig. 1 and the factor of 4 comes 

— * 

from forming the trace associated with the closed loop. We see that ps(k) does not depend 
upon k°, since we are making an on-shell approximation in the calculation of f(k — k'). In 
vacuum we have 

- iA(k) = -im° + (2G s i)ps(k) , (2.12) 

or 

A(k) = m° - 2G sPs {k) , (2.13) 

where ps(k) is real and negative. We remark when f(k — k') — 1, A(k) — > m and ps(k) — > 
(uu), so that we regain the usual result [6-8] 

m u = m° u -2G s (uu), (2.14) 

m d = m° d -2G s (dd) . (2.15) 

(Note that our G s is one-half of the G s defined in Ref. [7].) 

The evaluation of ps(k) proceeds by closing to contour in the complex k° plane. In the 
vacuum we find 

psCk) ^ 2N j^ . , (2 , 6) 



We define 



^ 3 ^k' 2 {l-B{k')) 2 + A 2 (k') 



E{k) = \f¥{\ - B{k)f + A 2 (k) , (2.17) 



and 



(k - k'Y = (E(k) - E{k')Y — (k — k'Y , (2.18) 

so that f(k — k') depends only upon \k'\ and the angle between k and k'. 

Integrals such as that in Eq. (2.16) require regularization. For our calculations we insert 
a factor exp[— k' 2 / a 2 ] with a = 0.60 GeV. (We have used the same Gaussian regulator in 
our calculations of meson spectra, where we have used a = 0.605 GeV [12-14]. However, 
those calculations included a model of confinement, so that we can not directly take over 



the parameters Gs, Gy and Gd used in those works.) We remark that an expression for the 
density-dependent condensate (uu) p may be obtained by using Eq. (2.16) with f(k — k') = 1. 

In the presence of matter it is useful to separate the propagator into two parts, one of 
which will give rise to the explicitly density-dependent terms. In this regard, it is useful to 
generalize Eq. (5.8) of Ref. [6]. We define 

E(k)^ - 7 • k{l - B{k)) + A(k) 



A (+) (£) 
A(-)(_jfe) = 



-E(k) 



2A(k) 

k(l-B(k))+A(k) 



2A(k) 



where E(k) = [k 2 (l - B{k)f + A 2 (k)] l/2 . Then 



S{k) = 



E(k) 



A(+)(jfe) AH(-jfe) 



k°-E+(k) k°-E-(k) 



+ 2m^-A( + \k) 9{kp - \k\) 5{k° - E + {k)) 
E(k) 

In the limit that A(k) — > m* and B(k) = C(k) = 0, we have, with 

E*(k) = 



S(k) 



m 



E*(k) 



A(+)(jfe) 



/c 2 + m* 2 , 
A(")(-jfe) 



2?r« 


m* 


E*(k) 




P- 


f m 


k 2 


— m* 2 



k°-E*(k) k° + E*(k) 
A( + \k)0(k F - \k\)5(k° -E*(k)) 



E*{k) 

which agrees with Eq. (5.8) of Ref. [6]. 
In the presence of matter, we have 



A(k,p)=m -2G s [pT(k)- P r(k)}, 



(2.19) 
(2.20) 



(2.21) 



(2.22) 
(2.23) 



(2.24) 



(2.25) 



where pg at (A;) is calculated in the same manner as p v g c (k), except that the upper limit of the 
integral over \k'\ is k F - Equation (2.21) is a generalization of a corresponding equation that 
may be found in Klevansky's review. (See Eqs. (5.18) of Ref. [6].) 

— # — # 

If we neglect B(k), Eqs. (2.11) and (2.16) provide a nonlinear equation for A{k) which 
may be solved by iteration. The results of such a calculation are reported in Table I where, 



for nuclear matter, we put kp = 0.268 GeV. In Figs. 2 and 3 we show values of the condensate 
(uu) and A(0, p) as a function of the density. These results may be usefully discussed in 
terms of the relation [15] 

H=H ( 1_ S + '") ■ (2 ' 26) 
Here is the pion-nucleon sigma term and p^ is the density of nucleons. If we put 

a N = 0.050 GeV, p N = (0.109 GeV) 3 , f n = 0.0942 GeV and m w = 0.138 GeV, we find a 38% 

reduction of the condensate at nuclear matter density, which agrees with our results given in 

Table I and Fig. 2. It is of interest to note that the linear dependence on the density implied 

by Eq. (2.26) appears to be valid up to about twice nuclear matter density. However, one 

may be concerned that, since we study quark matter rather than nuclear matter, Eq. (2.26) 

may not be appropriate. Consider, however, the relation 

H = Ho(l-^ + -) , (2.27) 

where a q is the quark sigma term and p q is the number density of the quarks, which we may 
put equal to 2>pN- Thus, if 3cr q = a at, we may use Eq. (2.26). The fact that a q ~ 15 MeV 
has been discussed by Vogl and Weise [7] . We have also discussed this matter in great detail 
in Ref. [16], where we calculated similar values of a q using the standard version of the NJL 
model. We conclude that the use of Eq. (2.26), or Eq. (2.27), with an appropriate value 
of cr q , is satisfactory. For example, if a q = cr N /3, as suggested in Ref. [7], the two relations 
imply the same density dependence of the condensate. 

In Fig. 3 we show A(0,p) which is the density-dependent mass parameter of the theory 
when B(k,p) = 0. We see that A(0,p) follows the trend seen in Fig. 2 for the density 
dependence of the condensate. 

III. LORENTZ- VECTOR TERMS OF THE QUARK SELF-ENERGY 

We now write £(£) = T, s (k) + Sy(A;) where 

= -7-£b(£)+7°C(£). 

For the calculation of C(k) we obtain the contribution from the last term in Eq 
with the result that 

C(k) = 2G vP v 2 {k) 
9 



(3.1) 
(2.21), 

(3.2) 
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FIG. 2: Values of the condensate (uu) are given as a function of 10 3 fefi. For nuclear matter 
10 3 fe| = 19.2 GeV 3 . [See Table I.] Here G s = 13.0 GeV~ 2 and B(k,p) is put equal to zero. 
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FIG. 3: Values of ^4(0, p) are given as a function of W 3 kp. [See Table I and the caption of Fig. 2.] 



with 



p V Ak) = 2N c j ^y 3 f 2 (k-k>) 



(3.3) 
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10 3 k 3 F 
( GeV 3 ) 


-(to) 1 / 3 
(GeV) 


-10 3 (uu} 
( GeV 3 ) 


A(0,p) 
(GeV) 





0.241 


14.0 


0.371 


10 


0.225 


11.4 


0.298 


19.2(nm) 


0.205 


8.67 


0.226 


25 


0.189 


6.75 


0.177 


30 


0.171 


5.00 


0.133 


40 


0.134 


2.46 


0.077 


50 


0.113 


1.44 


0.041 


60 


0.0990 


0.967 


0.030 


70 


0.0915 


0.898 


0.024 


80 


0.0853 


0.620 


0.021 


90 


0.0805 


0.522 


0.018 



TABLE I: Values of the condensate and A(0,p) are given for Gs = Gy = 13.0 GeV" 2 and 
m u = m d = 0-005 GeV. Here kp = 0.268 GeV for nuclear matter. We note the reduction of the 
condensate of 38% and a 40% reduction of A(0, p) at nuclear matter density [W 3 k F = 19.2 GeV 3 ]. 
Here a = 0.60 GeV is used in the Gaussian regulator exp[— k 2 /a 2 ]. 

An expression for B(k) may be found from the relation 

-i[-^kB(k)} = (-2G v i)(-l)i j^yS{k')f\k-k') (3.4) 

— * 

if we only keep the term proportional to 7 • k' in the expression for the quark propagator. 
In vacuum we may compare corresponding terms in Eq. (3.4) 

- 7 • kB(k) = -2G V / * »^f^ . (3.5) 

J (2tt) 4 [k> - E+(k>) + ie][k£ ~ E-(k') - it] 

Thus, 

B™ c {k) = 2G vP r(k) , (3.6) 
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with 

,r, /"^' \k%k-k')f 2 {k-k')\l-B(k')\ , „ N 

WpTW = 2N C / — — ' J >±±=J± LJJ . (3.7) 

W ^/£;' 2 [1 -S(ifc')] 2 + A 2 (A;') 

Here, fc and fc' are unit vectors. As in the calculation of A(k), Eq. (3.6) is generalized to 
read 

B(k) = 2G v [ P r(k)- P r(k)], (3.8) 

where pf at (k) is calculated using Eq. (3.7) with an upper limit on \k'\ of Uf- 

In Table II we present results of our calculation of the condensate, A(0,p), A(0,p)/[1 — 
B(0, p)\, C(0, p), B(0, p) and A(0, p) - A(0, 0). We also define 

U s (k,p) = A(k,p)-A(k,0), (3.9) 

where Us(k,p) is the density-dependent modification of A(k,0) in matter. 

It may be seen from the values given in Table II that there is a thirty percent reduction 
of the condensate at the density of nuclear matter, while the value of A(0) is reduced by 
thirty-nine percent. We would obtain a thirty percent reduction of the condensate if o^r = 39 
MeV. 

In Figs. 4 and 5 we exhibit values of A(k, p) and A(k, p)/[l—B(k, p)\ for various densities 
and in Fig. 6 we present values of C(k,p). Figure 7 shows the values of Us(k, pnm) and 
C(k,p NM )- 



IV. OFF-MASS-SHELL EFFECTS 

In our calculations we have neglected the k dependence of A, B and C. In this Section 
we provide some justification of that approximation. We can combine Eqs. (2.9) and (2.11) 
to yield a nonlinear equation for A(k°, k) in vacuum 

A{k\k) = -(2G s i)N c (-l) l^ ^f^iJ,^'} (4-1) 

J ( 27r ) - A(k'°, k') +ie 

Here we have neglected m° and B(k°, k). We see that k° dependence arises from f 2 (k — k') 
which is given by 

f\k-k')=exp[-(k-k') 2n /P}, (4.2) 
12 
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FIG. 4: Values of A(k,p) are given as a function of \k\ for various densities: a) 10 3 k F = [solid 
line]; b) W 3 k 3 F = 10.0 GeV 3 [dashed line]; c) W 3 k 3 F = 19.2 GeV 3 [dotted line]; d) W 3 k 3 F = 30.0 
GeV 3 [dot-dash line] and e) W 3 k F = 40.0 GeV 3 [short dash]. Here G s = 13.5 GeV" 2 , G v = 10.0 
GeV" 2 , m° = 0.005 GeV and a = 0.60 GeV. [See Table II.] 
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FIG. 5: The quantity A(k,p)/[1 — B(k,p)], which plays the role of a momentum- and density- 
dependent mass parameter, is shown. [See Table II and the caption of Fig. 4.] 
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io 3 4 

( GeV 3 ) 


-(uu) 1 ^ 
(GeV) 


-10 3 (nn) 
( GeV 3 ) 


A(0, P ) 
(GeV) 


A{0,p) 


U s (0,p) 
(GeV) 


C(0,p) 
(GeV) 


B(0,p) 


1—B(0 o) 

(GeV) 





0.2401 


13.85 


0.365 


0.347 








-0.0544 


10 


0.2279 


11.84 


0.292 


0.272 


-0.073 


0.0794 


-0.0736 


19.2 


0.2128 


9.64 


0.223 


0.203 


-0.142 


0.0988 


-0.1019 


30 


0.1755 


5.36 


0.114 


0.0942 


-0.252 


0.115 


-0.210 


40 


0.1441 


2.99 


0.0606 


0.0435 


-0.304 


0.126 


-0.394 


50 


0.1291 


2.15 


0.0432 


0.0279 


-0.322 


0.136 


-0.545 



TABLE II: Various values are given for the case G$ = 13.5 GeV 2 , Gy = 10.0 GeV 2 , = mP d = 
0.005 GeV and a = 0.60 GeV. Note a reduction of 30% for the condensate and 39% for A(0, p) at 
nuclear matter density, where k F = 0.268 GeV and lO 3 ^ = 19.2 GeV 3 . 




FIG. 6: Values of C(k,p) are shown for various densities: a) W 3 k 3 F = 10.0 GeV 3 [dashed line]; b) 
10 3 A; 3 = 19.2 GeV 3 [solid line]; c) 10 3 A; 3 = 30.0 GeV 3 [dotted line] and d) 10 3 *; 3 = 40.0 GeV 3 
[dash-dot line]. Here G s = 13.5 GeV -2 , Gy = 10.0 GeV~ 2 , m° = 0.005 GeV and a = 0.60 GeV. 
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FIG. 7: The values of U (k, pnm) = A(k, pNM) — A(k, 0) [solid line] and C(k, pnm) [dashed line] are 
shown. U (k, pnm) represents the density-dependent correction to the vacuum value of the scalar 
term of the quark self-energy. 

with 



(k-k'f = (k° - k' ) 2 - (k - k'y , 



(4.3) 



which differs from the on- mass-shell version given in Eq. (2.18). 

The solution of Eq. (4.1) for A(k°,k) is shown in Fig. 8 for various values of k°. It is 
seen that the dependence on k° is weak as was assumed in this work. 



V. THE NUCLEON SELF-ENERGY IN MATTER 

If one uses the Dirac equation to describe the interaction of a nucleon with a nucleus, 
or with nuclear matter, it is found that a strong scalar attraction is needed as well as a 
strong vector repulsion [17, 18]. The scalar field is of the order of -400 MeV and the vector 
field is about 300 MeV. It is of interest to see if the nucleon self-energy, Y, N = V s + 7°Vy, 
can be calculated in terms of the quark self-energy obtained in this work. To carry out 
this program we use a simple model of the nucleon in which a quark is coupled to a scalar 
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FIG. 8: Values of A(k°, k) are shown as a function of \k\ for various values of k°. 

diquark. The full complexity of the wave function, including vector diquarks and various 
relativistic effects, is discussed in Ref. [19]. 

We can calculate the nucleon self-energy in nuclear matter using a triangle diagram in 
which one of the lower two vertices of the triangle represents a vertex function for a zero- 
momentum nucleon to emit a quark of momentum k leaving a spectator (on-mass-shell) 
diquark of momentum —k. The other lower vertex represents the inverse process. At the 
upper vertex we insert the quark self-energy, = Us(k) + 7°C(fc) calculated in this work 
and integrate over k. We make use of Fig. 4 of Ref. [19] and parametrize the product of 
the vertex function and the quark Greens function by the quark-diquark wave function 

i,(k) = -^=e- p ' x2 u{k,s), (5.1) 



where u(k, s) is a spinor for a quark of momentum k and spin projection s [19]. 
The normalization for a single quark is obtained from the relation 

1 f d 3 k _ 2 p /X 2 ( k 2 \ _ 

Njwr e \ l + [E q ik )+mq Y)- 1 ' (5 ' 2) 
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where we put A = 0.18 GeV to correspond to the results of Ref. [19]. We may then relate 

— * 

the nucleon self-energy to the quark self-energy. For a nucleon of momentum P = 0, we 
have, with E q (k) = [P + m^] 1 / 2 , and m q = 0.364 GeV, 

*-*f&«^{ l + vd^)' (5 ' 3) 

and 

Vs = I fi^Usik) e'^ 2 ^ ( 1 S 1 . (5.4) 

The difference of sign in the brackets appearing in Eqs. (5.3) and (5.4) is due to the different 
behavior of the Dirac matrices, 1 and 7 , at the upper vertex of the triangle. 

Since the momentum content of the quark-diquark wave function is small [19], we expect 
that V s ^ 3E/s(0) and V v ~ 3C(0), so that V v ~ 296 MeV and V s ~ -426 MeV. A 
more careful evaluation of the integrals in Eqs. (5.3) and (5.4) yields Vy = 295 MeV and 
Vs = —392 MeV, which is in general accord with the values given in Refs. [17] and [18]. 



VI. DISCUSSION 



The behavior of matter at high density and low temperature has received a good deal 
of attention in the last few years [1-5]. A large part of the work in this area has made 
use of the NJL model. In our work we have attempted to modify the NJL model so that 
its predictions are in greater accord with QCD. As a first step, in Ref. [10] we introduced 
a momentum- dependent qq interaction which allowed us to reproduce the Euclidean-space 
behavior for the constituent quark mass obtained in lattice simulation of QCD [9]. In the 
present study we have considered the important vector interactions of an extended NJL 
model. Of particular interest is the behavior of the quark condensate at finite density. We 
find that the linear behavior in the density exhibited in Eq. (2.26) holds in our model up 
to about twice nuclear matter density. Also, we note that the use of a nonzero current 
quark mass is important for that result, since in the absence of an explicit chiral symmetry 
breaking term, the model exhibits a first-order phase transition at about 1.25 times nuclear 
matter density. 

Another point of interest are the results shown in Fig. 7. The quark self-energy is similar 
to that found in relativistic nuclear physics, with strong scalar attraction and strong vector 
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repulsion. The simple calculation reported in Section V suggests that the quark self-energy, 
when multiplied by 3, provides a satisfactory estimate of the nucleon self-energy which is in 
accord with results of relativistic nuclear physics [17, 18]. 

There is a body of work based upon the solutions of the Schwinger-Dyson and Bethe- 
Salpeter equations with a phenomenological form for the gluon propagator [20-23]. This 
body of work is reviewed in Ref. [24]. In contrast to the results of our work, it is found 
that the quark condensate increases with increasing chemical potential. The authors of 
Ref. [20] argue that the baryon density is zero up to the critical chemical potential, /i c , 
for deconfinement. They state that "This result is an expected consequence of confinement 
which entails that each additional quark must be locally paired with an antiquark thereby 
increasing the density of condensate pairs as fi is increased. For this reason, as long as 
ji < fi c , there is no excess of particles over antiparticles in the vacuum and hence the baryon 
number density remains zero." We note, however, that the baryon density usually considered 
in such calculations is due to the presence of nucleons, which have a finite baryon density. 
As noted earlier, in our work the baryon density is due to the presence of Fermi seas of up 
and down quarks, which serve to provide a model of the baryon density that would arise 
due to the presence of nucleons. We suggest that the conclusions presented in Ref. [20] do 
not refer to the situation of physical interest in which one studies the value of the quark 
condensate in nuclei or in nuclear matter. 
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